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Supersymmetry with heavy scalars is a model where at the LHC we have to rely on rate measure- 
ments to determine the parameters of the underlying new physics. For this example we show how 
to properly combine rate measurements with kinematic endpoints, taking into account statistical, 
systematic and theory uncertainties. Provided we observe a sizeable number of events the LHC 
should be able to determine many model parameters with small enough error bars to for example 
test unification patterns. 



I. SETUP 
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Supersymmetry as a prime candidate for new physics at the TeV scale should be discovered at the LHC even with 
moderate energy and luminosity The challenge for complex TeV-scale extensions of the Standard Model is to 
determine as many model parameters as possible at the TeV scale and extrapolate them to higher energy scales. 
This way, we can study the underlying structures and symmetries of an ultraviolet completion of our Standard Model 
■ — up to energy scales which might reach for example the scale of grand unification [EHU . 
Oh| Most studies which focus on understanding new physics at and above the TeV scale rely on a multitude of kinematic 
observables. In particular at the LHC kinematic measurements are the most powerful, because they can be extracted 
' in the presence of large QCD and top-pair backgrounds and are less prone to huge QCD corrections. Possible limits 
to such strategies we have seen in LHC studies of supersymmetry with light sleptons. There, the number of kinematic 
observables is drastically reduced and the remaining kinematic features do not determine the absolute new-physics 
mass scales well anymore The question then becomes how much information we can extract from fewer and less 
robust observables, including production rates of supersymmetric final states. Two aspects of such measurements mean 
O s \ additional complications: first of all, we do not actually measure a total signal cross section, but a matrix of production 
rates times branching ratios in the presence of backgrounds and possibly relevant kinematic cuts [lCf . Secondly, for 
such measurements the combination of experimental and theory uncertainties becomes the crucial stumbling block 
which determines if we can for example test gaugino mass unification at the LHC or not. This situation is somewhat 
similar to Higgs sector analyses at the LHC |ll| - [l3j . 



In supersymmetry with decoupled scalars (DSS), all scalar partners are decoupled from the relevant mass spectrum 
for the LHC 

m~ g = m q = m H ^ H ±^ = m$ > 10 4 GeV (1) 

This scalar mass scale might be very large [HI, [lj| or simply decoupled from LHC production above O(10 TeV) 
The observable spectrum consists of the usual Standard Model fields, the gluino g, the wino W, the bino B and the 
higgsino components H u c [- Omitting gauge-invariant kinetic terms and non-renormalizable operators, the Lagrangian 
of the low energy effective theory reads [lH, [13, [H| 
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where a a are the Pauli matrices (e = ia 2 ). It includes one light Higgs doublet, tuned to have a small mass m for 
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H = — cos (3eHj + sin /3H U . At the scale ms, the low-energy effective theory is matched to the full MSSM 
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The four parameters g are the Yukawa couplings of the neutralinos and charginos, which are modified with respect 
to their supersymmetric values. Since it will be impossible to observe these deviations at the LHC [l7| we do not 
include them in our parameter extraction. A t h are threshold corrections to the quartic Higgs coupling which change 
the tree-level Higgs mass m\ — 2Xv 
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with X t — A t — /i/tan/3 in the same range as the gaugino and higgsino masses. The weak-scale parameter tan/3 only 
appears in the boundary conditions and therefore is not a parameter of the low-energy effective theory. It is interpreted 
as the fine-tuned angle that rotates the two Higgs doublets into one heavy and one light mass eigenstate [l9j . 

New particles entering the renormalization group running at an intermediate scale ms {mz < ms < Mgut) 
contribute identically to the running of the three gauge couplings provided they compose complete representations 
of the unification group [2(|. All sfermions in the MSSM form complete SU(5) representations, so a possible gauge 
coupling unification scheme in the MSSM is unchanged by heavy scalars. Experimentally establishing such a pattern 
is one of the main long-term goals of the LHC. 

In our analysis we construct the universal gaugino masses at the GUT scale: M^Mqut) = J n \/p.- They are then 
evolved down to the scale mg based on one-loop renormalization group equations of the MSSM |2ll . [22j . The higgsino 
mass term /i is provided as an independent input parameter at the scale mz- Below ms we integrate out all scalars 
and run the modified renormalization group equations [15| to the desired scale. 



Keeping one of the supersymmetric Higgs bosons light requires a fine tuning which for this study we accept without 
offering an explanation. The Higgs mass matrix for the two Higgs scalars 
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has eigenvalues (m 2 H ) ± %/A 2 + b 2 in terms of (m%) = (m 2 H + m 2 Hd )/2 
the light Higgs mass to be of the order of the weak scale translates into 
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|/i| 2 and A = (rrijj — m 2 Hd )/2. Requiring 
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< (m 2 H ) < y/A 2 + b 2 . 
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The b term in the Lagrangian density breaks a Peccei-Quinn symmetry and can therefore be kept small, as opposed 
to A and (m 2 H ) which should both be of order mj. A light Higgs mass means that (m 2 H ) ranges around m 2 w . Hence, 
the fraction of the ((m 2 H ),A) space that satisfies Eq.© corresponds to 
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Heavy scalars leave the charging and neutralino spectrum untouched, so the lightest neutralino should still be 
a good dark matter agent [l5l |23j |. The measured density of dark matter then imposes a constraint on /i, which 
unlike in mSUGRA toy models is not determined by electroweak symmetry breaking. The WMAP measurement 
JJdm'' 2 = O-Hlio.008 ^3] can be reproduced in different parameter regions, which will be represented by our choice 
of reference parameter points: 

- the 'mixed region' with Mi w \i and a mixed higgsino-gaugino LSP. Here, XiXi annihilation is enhanced for 
gauge and/or Higgs bosons or top quarks in the final state. 
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10 TeV 


m 1/2 (M GVT ) 


132 GeV 297 GeV 


fj,(m z ) 


290 GeV 200 GeV 


tan /3(ms) 


30 


A t {m s ) 






TABLE I: Left: two parameter points chosen for this LHC s~ 
SuSpect All masses in GeV. 





DSS1 


DSS2 
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129 
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60 
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185 


X° 3 


296 


206 


X°4 


310 


317 


xt 


117 


175 


xi 


313 


317 



'. Right: relevant mass spectrum from a modified version of 



- the 'pure higgsino' and 'pure wino' regions where the LSP is almost mass degenerate with the Xi an( i the x.2- 
This leads to an enhanced co-annihilation. This region generally requires an LSP heavier than 1 TeV. 

- the 'Higgs pole' region in which the LSP is rather light, to^o rs \m,h and the annihilation proceeds via a resonant 
light Higgs. 



If gluinos are lighter than squarks, they will mainly decay through virtual squark exchange into quarks and 
charginos/neutralinos [25| . For very heavy squarks quantum corrections to the gluino decay processes can be sig- 
nificant, because they are enhanced by the large logarithm m?,/m q 1261) . If the scalar mass scale is larger than 
10 6 GeV, the gluino becomes sufficiently stable to form i?-hadrons [27| which can be analyzed at the LHC without 
major difficulties [l7], HI]- Those are subject to cosmological constraints if they affect nucleosynthesis. A gluino with 
TeV-scale mass must have a lifetime shorter than 100 seconds to avoid altering the abundances of deuterium and 
lithium-6. This sets an upper limit of mg < 10 9 GeV [29| . Therefore, in this study we will focus on comparably 
short-lived gluinos. 

To quantitatively study supersymmetry with heavy scalars at the LHC we define the two parameter points shown 
in Table HI They are in agreement with constraints from both dark matter observations and collider searches. A 
scalar mass scale of 10 TeV is the lowest value which still qualifies as 'decoupled' on LHC energy scales. It limits the 
amount of fine-tuning in the Higgs sector to one part in 10 4 and gives us a short enough gluino life time to avoid 
non-standard phenomenology or undesirable cosmological effects. 

The most important model parameters are \i and mi/2- They define the mass spectrum as well as the field content 
of the neutralinos and charginos and most notably the LSP. Our two choices span most of the LHC-relevant parameter 
space allowed by LEP and WMAP. DSS1 lies in the Higgs pole region, so the LSP is mostly a bino with a non-vanishing 
higgsino component. DSS2 lies in the 'mixed region'. The light LSP leads to an invisible Higgs branching ratio around 
1%. In DSS1, the gluino, the x\ and xf are fairly light and have correspondingly large LHC cross sections, as shown 
in the Appendix. For the heavier DSS2 spectrum this is no longer the case. 

The large value tan (3 = 30 avoids LEP limits on the light Higgs mass. Its only impact on the low-energy theory is 
on m,h- The effect of A t on rrih is suppressed by at least l/m|, so we can as well set A t = 0. 



II. LHC OBSERVABLES 



Throughout our analysis we use leading-order and, where available, next-to-leading-order Monte Carlo generators 
for signal and background processes. The LO-generated samples are normalized to NLO cross sections using K- 
factors. For SUSY and VV events we rely on HERWIG [3(| including initial/final state radiation, spin correlations in 
the decay of heavy states and angular correlations between jets. The V + jets events we obtain from ALPGEN 
including matrix element and parton shower matching [321 . T he NLO normalization of these rates is given by MCFM [33] 
or PR0SPIN02 [34]. Top pairs we simulate with MCONL0 [35j . including its NLO normalization. The pure QCD jets 
background we expect to be heavily suppressed by cutting on, e.g. missing energy. Therefore, we leave it at a leading- 
order PYTHIA simulation [36j . Detector and reconstruction effects we account for with a standard general purpose 
detector simulation as described in Ref. 181. 



Looking at the DSS signatures we rely on standard observables. They are based on the reconstruction of isolated 
jets with R = 0.4 and leptons with px > 20 GeV. From the missing energy measurement we can determine the 
effective mass M c h as a measure of the total activity in the event. It includes the four hardest jets and all identified 
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after cuts 


add'l M cS > 






800 GeV 1 TeV 


DSS1 


12631 


4701 


DSS2 


145 


112 97 


tt 


5161 


274 62 


QCD jets 


848 


15 


VK+jets 


769 


195 


Z+jets 


422 


128 


WZ 


20 


4 




1 9 


Q 
O 


ZZ 


2 


2 


total SM 


7234 


621 62 


DSS1 significance 




18 


DSS2 significance 




0.8 6 



TABLE II: Number of events remaining for each process after the cuts listed in the text. They are normalized to NLO cross 
sections for an energy of 14 TeV and an integrated luminosity of 1 fb _1 . 

leptons. For SUSY events, M e a scales with the mass of the heavy particles produced and can be used to quantify 
the mass scale of SUSY events [33, EH]- The transverse sphericity is defined as St = 2A 2 /(Ai + X 2 ) where A, are the 
eigenvalues of the 2x2 sphericity tensor = ~}2 k PkiPkj- This tensor we compute using all jets above pt = 20 GeV 
and all leptons. SUSY events tend to be relatively spherical (St ~ 1) since the initial heavy particles are usually 
produced approximately at rest and their cascade decays emit particles in all directions. In the DSS1 parameter point 
the sparticles are fairly light, resulting in a uniform distribution between St = 0.1 and 0.6. For QCD jets or U+jets 
events St peaks at zero. A cut on St does not reduce the tt background as the distribution is very similar to the 
signal. 

The most frequent final states occurring for the DSS1/DSS2 parameter point can be classified by the number of 
leptons: U (70%/70%), U (20%/23%), 21 (5%/5%), and 3£(3%/l%). The number of decay jets can vary from zero 
in the case of leptonically decaying xfx.2 pairs to eight in the case of gluino pairs. This does not include jets from 
the underlying event or initial and final state radiation which have to be carefully studied in addition 39]. The 
most common channel for both points is the zero-lepton channel gg — > xfxt + 4 jets — ► 8 jets + We use it to 
estimate the discovery potential using standard cuts [38;]: One hard jet with pt > 100 GeV and three over 50 GeV, no 
electrons or muons, an effective mass M c g > 800 GeV, missing transverse energy $t > 100 GeV with fx > 0.2 x M e s, 
sphericity St > 0.2 and finally a geometric separation between jets and the missing energy vector of A</> > 0.2 for the 
three leading jets. Both for jet mis-measurement and b decays the fx vector will be close the direction of one jet, so 
this cut reduces fake missing energy from QCD. Note that for the various LHC observables discussed below this basic 
set of cuts will be modified. 

Table [TT1 shows the number of events remaining after cuts. After all but the M e g cut tt is the dominant background, 
but there are also significant contributions from U+jets. Finally, M e g > 800 GeV reduces the background to below 
the level of the signal for DSS1. 

The systematic uncertainties on the number of background events for 1 fb _1 we take to be 50% for QCD multi-jets 
and 20% for tt and U+jets, WW, WZ and ZZ. This corresponds to a combination of data-driven and Monte Carlo 
methods (38j . 

A discovery of new physics can then be claimed if the number of observed events exceeds 25 and the significance 
is larger than five. The significance of the observation of DSS1 with an integrated luminosity of 1 fb _1 is 18, so 
this parameter point will be discovered at the LHC within one year of data taking at low luminosity. For DSS2, the 
significance is low and, for this set of cuts, does not increase with statistics. However, if we require M c g > 1 TeV the 
significance for 1 fb _1 increases to 6. 

Higgs mass 

With decoupled scalars the light Higgs scalar is essentially equivalent to its Standard Model counter part at 129 GeV. 
Its mass depends on Ms and on tan /3. The Ms dependence arises from the running of A from ms to the weak scale, 
while tan j3 appears as cos 2 2/3 in the matching. This can impact the numerical value of rrih by up to 20 GeV and 
should — like usually in supersymmetric Higgs studies — allow us to determine tan/3 at the LHC. 
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theory 


fit value 


statistical error 


DSSl 






55.1 


55.2 


±0.6 / 10 ftr 1 


DSS2 




-m x o 


60.7 


60.2 


±2 / 100 fb" 1 




x 3 




81.9 


79.0 


±3 / 100 fb" 1 



TABLE III: Results of the fit to the invariant mass distribution. All values given in GeV. 



The total next- to- leading order production cross section is 39 pb. It is computed by HIGLU for gluon fusion [40| . 
VV2HF for weak boson fusion and V2HV for the production in association with a vector boson. This number includes 
NLO QCD corrections, and for the first one also the NLO electroweak contributions [4l|. The dominant decays are 
into bb (53%), WW* (29%), rr (5%), and ZZ (4%) [H|43j]. In addition, the branching ratio into photons reaches its 
maximum of 0.2%, allowing for a precise mass measurement. 

Higgs production through supersymmetric cascades only occurs in x§ 4 an d xf decays. The corresponding pro- 
duction rates are small, so their contribution to the total Higgs production is of the order of 100 fb, i.e. negligible 
compared to SM channels. 

Systematic uncertainties on the Higgs mass measurement arise from the electromagnetic energy scale. The calibra- 
tion of the photon energy scale will be derived from Z — ¥ ee events and Z — > fj.fX'y events, with an expected accuracy 
of 0.1%. The statistical uncertainty for an integrated luminosity 100 fb -1 should also range around 0.1% [381 . T he 
theory uncertainty due to higher-order corrections to rrih should not exceed the very conservative limit of 4% [44] . 



Di-lepton endpovnt 

Kinematic endpoints are usually the main ingredients to supersymmetric parameter analyses, due to their small 
experimental and theory errors [45144 7| . Perfect triangular di-lepton edges in cascade decays occur in successive two- 
body decays, like X2 ^ Xi^- Such a measurement directly constrains the gaugino mass parameters Mi and M2 
as well as the higgsino mass parameter \i. In contrast, decays via on-shell Z bosons only give us the Z mass peak 
with little information on the supersymmetric masses in the decay. For our parameter choices, x§ (DSSl) or both 
X3 4 (DSS2) as well as the charginos in DSS2 decay through such an on-shell gauge boson [43|.|48T|. 

However, we can also apply endpoint techniques to three-body decays: x\ ~^ Xi^ m both points and X3 —> Xi^ m 
DSS2 lead to two leptons with opposite signs and same flavor (OSSF), as listed in the Appendix. On the production 
side, the X1X2 channel has little background due to the small number of jets in the final state. To increase the total 
rate we also include X2 production from gluino decays. The total available cross section leading to this decay becomes 
ct(x° X°M) ~ 3.5 pb and 93 fb in DSSl or DSS2 and a{ X " -> x\ll) ~ 75 fb in DSS2. 

In addition to the staggered jet cuts (pr.j > 100, 50, ... GeV) we now require at least two OSSF electrons or muons 
withpx > 20 GeV, \rj\ < 2.5, and m« < m x o — m x o +10 GeV. Since the true value of the endpoint is a priory unknown, 
this choice implies that the edge has already been observed. To remove combinatorial as well as top backgrounds we 
apply flavor subtraction. Many backgrounds cancel in the combination 

N{e+ p E ~ ] +PN( i i + ^)~N(e ± ^) (8) 

where j3 is an efficiency correction factor equal to the ratio of the electron and muon reconstruction efficiencies. 

In the DSS2 parameter point, the lowest endpoint corresponds to the mass splitting m % o while the second one 
corresponds to m^o i . We fit its distribution with a superposition of three components, two modeling the decay 
kinematics plus a Breit-Wigner Z line shape. Table Mil compares the results with the theoretical values. The 
extracted values are in good agreement with the input values for DSSl and in reasonable agreement for DSS2. The 
statistical errors on the mass differences we can extract from the fit to the reference function. Systematic uncertainties 
are dominated by the lepton energy scale (0.1%), while theory errors due to unknown higher-order contributions are 
expected to range around a percent. 

There might occur doubts if the second edge giving m x o i is actually visible, so we check that the curve between 
the m x o edge and the onset of the Z peak indeed lies 5a above the background-only prediction. 

Assigning these measured values to sparticle mass differences necessitates a few assumptions. In DSSl, lepton pairs 
are quite frequent with respect to the overall SUSY production. This suggests that the neutralino triggering this 
decay is somewhat light. In addition, a decay through a 296 GeV x° is rather unlikely, and in such a case additional 
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DSS1 


a 


BR 




DSS2 


a 


BR 


<t(3£) 


direct 


11.7 pb 


1.54% 


180 fb 


direct 


1390 fb 


1.54% 


21.4 fb 


via gg 


10.4 pb 




160 fb 


via gg 


166 fb 




2.6 fb 




total 


340 fb 


total 


24 fb 



TABLE IV: Cross sections contributing to the tri-lepton signal for the LHC running at 14 TeV. 

structure would be seen. In DSS2, with two endpoints and a Z peak, the interpretation is more complicated. In 
addition to the assumption that the endpoints arise from X23 decays we have to assume that x° decays preferably to 
Xi- Otherwise, the largest endpoint could correspond to m x o and the Z peak to the decay ~^ Xi%- 

Di-jet endpoint 

Unfortunately, the technique described above is only applicable to decays involving b jets, but not light-flavor jets. 
In DSS1, 1.7% of gluinos decay to the LSP with two bottom quarks. We select these events by requiring four jets 
with Et > 50 GeV, no leptons and > 100 GeV. For two 6-tagged jets we compute the invariant mass rribb- The 
background is dominated by tt events, as well as combinations due to decays other than g — > Xi^b- The fit output 
from rribb we compare to the theoretical values: for 10 fb -1 the fit value of 380.6 GeV corresponds to the input of 
383.0 GeV within the statistical error of 5.2 GeV. This measurement is the basis for the determination of M3 from 
the gluino mass (46|. 

Tri-lepton cross section 

In contrast to the usual and more optimistic scenarios @, [46|, |4j| , decoupling all scalars at the LHC implies that we 
will not have enough kinematic information to extract masses and model parameters of the underlying new-physics 
model. Therefore, we need to rely on cross section measurements, in spite of their larger experimental and theory 
errors. The main purpose of this analysis is to show how such rate measurements can indeed be used as input to new 
physics measurements. 

As usually, signatures involving leptons have lower LHC backgrounds and increase the precision of the measurement. 
For heavy supersymmetric scalars charginos and neutralinos will give different final states with numerous isolated 
leptons. The tri-lepton final state allows for background rejection by requiring two OSSF leptons. It arises from 
chargino and neutralino production with subsequent decays X2 ~~ ^ Xi^ an d xf -» xliiy- In Table [TV] we show the 
composition of the tri-lepton signal in DSS1. Pairs of X1X2 are produced directly as well as in gluino decays. 

We select these events requiring at least one OSSF pair and exactly three leptons. In case of direct production, 
two LSPs will be emitted essentially back-to-back, hence canceling the missing transverse energy, so we lower our 
cut to fx > 50 GeV. An optional jet veto above pr ~ 20 GeV (dependent on detailed QCD studies) can be applied 
in order to select events from direct X\X\ production rather than g pair decays. However, its effect on signal and 
background rates is hard to predict, so the results should be taken with a grain of salt. To reject Z decays we veto 
mu = 81. 2. ..102. 2 GeV. For our signal, we expect m« < 56 GeV due to the x° 1 mass splitting. 

In Table IVl we present the number of events after cuts. Only tt and WW/WZ/ZZ are significant, where the latter 
is already partly removed by the mu cut while the former can be removed by a jet veto. 

An as precise as possible extraction of the number of tri-lepton signal events relies on our knowledge of the back- 
grounds and a complete understanding of detector effects, luminosity, parton distributions, and finally cut efficiencies. 
The systematic uncertainty on er(SUSY — > 3£) is certainly bounded from below by the knowledge of the luminosity 
C. i.e. of the order of 5% [50(. In order to take into account additional systematic errors we consider cases of 5%, 
10% and 20%. The theory uncertainty due to QCD effects we estimate to be of the order of 12% [34] ]. 

Gluino pair cross section 

The gluino pair channel constitutes a large fraction of the signal events in both parameter points (77% in DSS1 
and 22% in DSS2). Different strategies should allow us to select only gg events. However, they will be very model 
dependent, as they require at least a guess of the gluino decays branching fractions. We can for instance take 
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advantage of the very short zero-lepton cascade gg — > x\jQqqqq to remove SUSY background from chargino and 
neutralino channels. Again, the systematic uncertainty is bounded from below by the knowledge of the luminosity, so 
we consider systematic errors of 5%, 10% and 20%. Theory uncertainties from higher order contributions can be 30%, 
estimated from a scale variation by a factor 1/4... 4 at next-to- leading order [34J. It is meant to be conservative to 
also accommodate additional sources of uncertainties, like parton densities and the strong coupling. It can be reduced 
once we systematically include higher-order QCD corrections to the production of heavy particles |5l|. This choice of 
theory errors will give us a conservative estimate if including rate information at the LHC should become a part of 
supersymmetric parameter analyses. 



On-shell vs off-shell Z bosons 



If the mass difference between two neutralinos or two charginos is larger than mz, the invariant mass distribution 
of its decay products will exhibit a very sharp peak. In contrast, if the mass difference between the two sparticles is 
too small, the mu distribution will show a triangular shape with a sharp endpoint below mz- This effect can provide 
valuable information about the neutralino and chargino masses as well as their couplings to the Z boson. 

The measurement of the er(SUSY — > Z) rate requires a good knowledge of the luminosity, lepton efficiencies and 
background rates. Some of these source of systematic uncertainties cancel from the ratio 

Rz _ N (m u > endpoint) 
N [mu < endpoint) 

For our two reference points we hnd Rz < 0.004 ~ (DSS1) and 0.196 (DSS2), respectively. Systematic errors arise if 
the pt spectra of leptons from the Z peak and from the triangular shape are different or if the identification efficiencies 
for these two spectra are different. However, electron identification efficiencies are reasonably flat for pt > 25 GeV. 
Leptons from the Z peak will have transverse momenta around 45 GeV, while those under the triangle curve we 
cut to pt > 20 GeV. We therefore include an overall 1% systematic error to mainly account for lepton identification 
uncertainties. The theory uncertainty due to the prediction of branching ratios should be of the order of 1% jH, l48j. 



III. PARAMETER DETERMINATION 



Once new physics (e.g. supersymmetry) will be discovered at the LHC, we have to turn our focus towards under- 
standing the corresponding signatures. On the one hand, we will have to test different types of TeV-scale models with 
the available data, while on the other hand we have to consistently determine the parameters of the underlying theory 
for each of these model hypotheses. This might well include combining LHC observables with other measurements 
such as the relic density of dark matter, the magnetic moment of the muon, or flavor physics. Note that a consistent 
approach does not allow for the replacement of some measurements by top-down predictions in someone's favorite 
model. Instead, we need to see how far we can get, for example with the limited set of observables described in the 
last section. 

SFitter 0, H2 is designed to map up to 20-dimensional highly complex parameter spaces onto a large set of 
observables of varying quality, which can be highly correlated. It can be used to estimate the reach in terms of a 
given model for any experiment, but also to realize a proper bottom-up approach to determine the parameters of a 
fundamental theory. 

The determination of the parameters then proceeds in two steps. First, we maximize the exclusive log-likelihood 
using a weighted Markov chain [2] to identify the best-fitting point in parameter space. The starting point of this 
Markov chain is arbitrary, and we repeat the search several times to ensure our procedure converges well. This 





after cuts 


add'l jet veto 


DSS1 


681 


43 


DSS2 


87 


4 


tt 


1,106 


59 


QCD 








W/Z + jets 


14 





WW/WZ/ZZ 


235 


73 



TABLE V: Number of tri-lepton events remaining after the cuts discussed in the text (assuming 10 fb 1 ). 



8 







uncertainties (%) 




observables 


stat. 


systematic 


th. 










~ i 
value 


source 




DSS1 


m h 


129 GeV 


0.1 


0.1 


energy scale 


4 




X 2 Xi 

m-g - m x o 


oo. z Lrev 


1 


U.l 


energy scale 


1 




382.8 GeV 


1.5 


1 


energy scale 


1 




a\il) 


340 fb 


2 


> 5 


luminosity 


12 




Rz 


< 0.004 


0.01 


1 


lepton id. 


1 




a (gg) 


DZ.O pU 


0.1 





luminosity 


30 






129 GeV 


U. 1 


0.1 


energy scale 


A 

4 




m„o — m„o 

X-2 Xi 


60.7 GeV 


3.3 


0.1 


energy scale 


1 




<r(a«) 


81.9 GeV 


3.7 


0.1 


energy scale 


1 




23 fb 


14 


> 5 


luminosity 


12 






0.57 


0.7 


1 


lepton id. 


1 




a{99) 


1067 fb 


3 


> 5 
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TABLE VI: Left: summary of available collider observables in DSS1 and DSS2. Statistical errors are quoted for an integrated 
luminosity of 100 fb -1 . Right: absolute and relative errors on the determination of the underlying model parameters for three 
fitting strategies described in the text. 



minimum then serves as starting point for a MINOS hill-climbing minimization to improve the resolution and to 
estimate the errors. 

In models with decoupled scalars we can use the different LHC observables discussed above to determine the 
parameters of the model. Table IVII summarizes them along with their expected uncertainties for both parameter 
points. For the mass differences we use the result of detailed experimental analyses, while for the rate-related 
observables we rely on the theoretical central value, lacking the complete experimental analysis. 

By definition, no information on the squark and slepton sector is available except for its explicit absence. Conse- 
quently, we fix ms and A t to large (nominal) values. The three gaugino mass parameters we fit independently, to 
allow for a bottom-up experimental test of gaugino mass unification. Technically, we know that for scalar masses a 
consistent bottom-up approach does not reproduce the usual top-down results, which means we would have to evolve 
all parameters strictly from the weak scale to the high scale (jf . For gauginos the differences between the two methods 
are not as large, so for illustration purposes we use a top-down running for the renormalization group equations. 

Four scenarios illustrate well the precision we can reach on the determination of supersymmetric model parameters 
with decoupled scalars. We use about 1000 toy experiments for each scenario and each benchmark point. The toy 
experiments are generated by smearing the observables according to the expected experimental and/or theoretical 
errors, depending on the scenarios. Correlations among the measurements are taken into account separately for the 
energy scale of leptons, jets and the luminosity measurements. For each toy experiment the best-fit parameter set is 
determined. From the distribution of the best-fit parameter we read off the error as the RMS (Root-Mean-Square) of 
the distribution. Table IVTl displays the resulting uncertainties for the four scenarios: 

- EXPerimental errors. In order to evaluate the impact of the pure experimental uncertainties on the determination 
of the DSS parameters for a given luminosity of 100 fb -1 , we take into account the statistical and systematic 
errors, but we do not include theory errors. 

We see that with roughly a 50% error on tan /3 we can hardly determine this parameter in both reference points. 
Because all heavy Higgs bosons are decoupled our only leverage is the light Higgs mass which depends on several 
parameters, including mg in the (s)top sector. The better way to study the Higgs sector, possibly including 
tan/3, would be a dedicated Higgs analysis at the LHC [HI, [HI. 

The second Higgs-sector parameter fx is determined to better than 1% in DSS1 and 4% in DSS2. In the former, 
this is due to large production rates for neutralinos and charginos. In the latter, the higgsino fractions are 
well spread over all neutralinos, so neutralino mass differences include this information. Thanks to very small 
statistical uncertainties on all mass differences in DSS1, also the gaugino mass parameters are determined within 
2%. This is not the case in DSS2, implying a deterioration to ~ 10%. As discussed above, the observability 
of the second neutralino edge is arguable. We perform the DSS2 fit with and without this observable and find 
identical results. 

- SYSTematic errors. Here we assume that a very large number of events has been gathered at the LHC (> 
300 fb -1 ) and that statistical uncertainties are negligible. By that time, we might assume that theoretical 
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predictions of nih and the different cross sections will have rendered the theory error negligible as well. This 
idealized scenario is useful to measure the impact of systematic errors and the ultimate precision of the parameter 
determination. 

In Table fVTl we observe that with only systematic errors all parameters are determined within 1%, except for 
tan/3 which still suffers from an invisible extended Higgs sector. Systematic errors on cross section measurements 
we vary between of 5 and 20%. In DSSI, this hardly affects the parameter determination as all model parameters 
are already well constrained by mass measurements. In DSS2, a 20% systematic error on a{gg) doubles the 
uncertainty on M3 as compared to the case where we are dominated by luminosity measurement (5% systematic 
error). In contrast, a 20% systematic error on cr(3£) again does not affect the parameter determination, because 
the weak masses are constrained by kinematics measurements. 

Compared to the first case including statistical and systematic uncertainties we observe a gain of at least a 
factor 4 in the errors for DSS2. This shows that for the integrated luminosity of 100 fb _1 assumed in the 
scenario EXP the statistical error still dominates. 

- THeory errors. Again, we assume very large statistics but include theory errors, thus taking into account 
theoretical and systematic errors. This scenario gives us a flavor of what is achievable after at least five years 
of operation of the LHC at full luminosity. 

With a 4% theory uncertainty on mu, tan j3 is now practically undetermined. The interplay of neutralino mass 
splittings, Rz and the tri-lepton cross section provides enough constraints on the neutralino and chargino sector 
to allow for a determination of Mi, M2 and \x to better than 4%. In DSSI, M3 is very well determined by the 
mass splitting m~ g — m x o . However, in DSS2 the only available handle on M3 comes from the gluino pair cross 
section which suffers from a large theoretical uncertainty which we set conservatively to 30%. As the production 
cross section strongly decreases as function of M3, the distribution of the toy experiments resulting from the 
symmetric theoretical error is highly asymmetric with a long tail to large values of M3 as expected. In spite of 
the large theoretical error on the cross section prediction, a M3 measurement significantly better than at the 
10% level is feasible. 

- FULL errors. In this scenario we combine the experimental errors (EXP) with the theoretical errors (TH) on the 
observables for 100 fb^ 1 . In DSSI, the errors on the gaugino masses are slightly larger than the errors for the 
TH scenario. In DSS2 the errors on Mi and M2 are essentially the same as the EXP errors, the corresponding 
mass difference measurements are dominated by the statistical error. For M3 the experimental errors are fairly 
small compared to the theoretical error. Nevertheless, the determination remains at the level of several percent 
due to the steep descent of the cross sections as a function of the gluino mass. Including further theoretical 
developments is necessary to increase the precision of the parameter determination by reducing the theoretical 
error 15 ill. 



IV. OUTLOOK 

Supcrsymmctry with heavy scalars is a variation of the MSSM which more naturally accommodates for example 
flavor constraints, at the expense of the solution of the hierarchy problem. For the LHC, it is irrelevant as which 
energy scale the scalars reside, as long as they are at least of the order of 10 4 GeV. 

Such a model is a serious challenge to any kind of supersymmetric parameter analysis, because it severely reduces 
the number of LHC observables, in particular from cascade decay kinematics. Instead, we need to rely for example 
on rate measurements, including their complex systematic and theory error structure. Of our two parameter points 
the first one should be discoverable within a year of data-taking at the LHC, due to very large SUSY cross sections. 
The second point has lower rates but should be discovered within a few years. For both scenarios we establish a set 
of observables, including statistical, systematic, and theory uncertainties. 

For a model with heavy scalars we have shown that a global fit of the model parameters to the experimental 
observables is still able to determine the correct central values and corresponding errors for all parameters. The 
weakly interacting sector (Mi, M2 and /i) can be fairly well measured at the LHC, with conservative accuracies of 
the order of a few percent after including all error sources. Most notably, this includes all uncertainties related to 
rate measurements at the LHC. The Higgs sector suffers from the fact that we will only have one observable at hand, 
namely the light Higgs mass. For such a situation we will have to resort to a dedicated Higgs sector analysis at the 

lhc dim. 

Even for optimistic LHC luminosities and energies the impact of systematic uncertainties is likely not dominant, 
even though an estimate of these errors prior to a full-fledged analysis on real data should be taken with a grain of 
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TABLE VII: Next-to-leading order cross sections for direct production of sparticles at the LHC for DSS1 and DSS2 as computed 
by Prospino2 The individual rate add to 81.6 pb (DSS1) and 4.3 pb (DSS2). 

salt. For all theory uncertainties we have relied on particularly conservative estimates, which means that the positive 
outcome of our study is generally dependable. 
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Appendix A: Production rates and branching ratios 

In contrast to most supersymmetric parameter analyses for the parameter point SPSla [491 ] . the analysis presented 
in this paper heavily relies on rate information. Note that 'rate' refers to cross sections times branching ratios for 
signal plus background. Therefore, we list all cross sections as well as all branching ratios computed at NLO in this 
appendix. 

Table \VU\ lists the next-to-leading order cross sections for direct production of supersymmetric particles at the LHC 
for DSS1 and DSS2 as computed by Prospino2 [34| . The combined production rate for strongly interacting particles 
suffers in models with decoupled scalars most of all because usually the quark-gluon initiated associated squark-gluino 
channel dominates the supersymmetric LHC samples. Gluino pair production is still more than a factor two larger 
than light- flavor squark-antisquark production, so we can still expect sizeable SUSY samples at the LHC. For the 
DSS1 parameter point the light gluino indeed yields a relatively large cross section for the production of gluino pairs. 
In DSS2, gluinos are fairly heavy and the gluino production rate is low. 

The second largest contribution to SUSY production is due to the associated production of charginos and neutralinos 
with 12 pb and 2 pb in DSS1 and DSS2, respectively. In DSS1, this channel is completely dominated by xtx® 
production while in DSS2, equivalent contributions arise from xtxl 2 3 an d xtxl- This is because the neutralino 
mass splittings are smaller and the gaugino/higgsino content is more degenerate in the heavier DSS2 parameter point. 
Its lower value of /1 increases xtx.3 production but fails to keep up with the very light Xi 2 m DSS1. The only other 
significant contribution to SUSY production in the light DSS1 point comes from xt pahs while in DSS2 there are 
contributions from neutralino pairs. 

Table IVIIII lists the higher-order branching fractions of the supersymmetric particles involved. Most gluino decays 
contain three particles in the final state, proceeding through a very virtual squark. Due to its wino-like nature, the 
X2 decays like a Z boson plus missing energy. The larger higgsino component in the DSS2 parameter point brings in 
decays to charginos. Again, the light chargino decays like a W boson plus missing energy. 
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TABLE VIII: Branching fractions for both scenarios computed by SDECAY [4a.l48l|. Values rounded to the full percentage. 



Appendix B: Error treatment 



As discussed in this paper, the focus on rate measurements in the case of heavy scalars forces us to generally treat 
the theory uncertainties on rate observables as possibly the dominant error at the LHC. At the same time, some of the 
production rates are small, so our Markov chain might run into parameter regions with small event numbers which 
have to be treated using Poisson instead of Gaussian statistics. 

Like in all SFitter analyses 0, [l2| we follow the Rfit scheme [52[ to combine Gaussian experimental and flat theory 
errors. This scheme interprets theory errors as a lack of knowledge on a parameter. As long as the deviation between 
theory and experiment is within the theory error, this must not have any influence on the total likelihood. Once 
the difference becomes larger the (perturbative) theory is simply ruled out instead of just very, very unlikely. This 
means that we cannot simply convolute some kind of theory likelihood distribution with a Gaussian experimental 
error. Even assuming a flat theory likelihood curve would give the difference of two one-sided error functions, i.e. a 
peaked likelihood. Instead, the combined log-likelihood which can also be derived using a profile likelihood ansatz is 
given by 



for \di-di\< of 

for \di -di\> a\ 



21ogL = X 2 = { m-dil-ofY th ■ ( B1 ) 



For large enough event numbers the experimental error is a combination of three different sources, all Gaussian and 
summed in quadrature. The statistical error is uncorrelated between different measurements. A first systematic error 
originates from the lepton energy scale and the second from the hadronic energy scale. They are treated separately. 
Each is taken as 99% correlated between different observables. 

For smaller signal numbers systematic and statistical uncertainties are incorporated by convoluting a Poisson 
probability with a Gaussian probability, where the number of background events Nf, is the mean and 6b (systematic 
uncertainties) the standard deviation [53|. The probability that the background fluctuates at least to the observed 
number of events A^ b s = Signal + M> is 

p = A d6 Gauss(AMb) V — (B2) 
■'• ^ J] 

where A normalizes the integral. Then the significance of the signal reads Z n — -\/2erfo 1 (l — 2p). If N \, s is very large 



compared to Nb, this significance is approximated by 54 1 




Z Q = —==Z== I W N ohs + | - J-^Vfe + ||r I • ( B3 ) 



VT+Sj/Nb \ V '8 y 8N b 
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